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ABSTRACT 

We argue that gravitational instability in the outer parts of collapsar disks may lead to fragmentation near the 
radius where helium photodisintegrates, because of the strong cooling provided by this process. This physics 
sets clear physical scales for the fragmentation conditions and the properties of gravitationally bound clumps. 
Collapse of a fragment proceeds until the neutrons become degenerate; a neutron star of mass wO.I-IMq 
may result. We find that tidal disruption of a fragment and accretion by the central black hole are too rapid to 
account for the durations of observed X-ray flares from long gamma-ray bursts. Prior to disruption, migration 
of the fragment is driven by gravitational radiation and disk viscosity, which act together to produce a unique 
gravitational-wave signature. Advanced LIGO may be able to detect such sources within wlOOMpc. 
Subject headings: accretion disks — black hole physics — gamma rays: bursts — gravitational waves 



1. INTRODUCTION 

Observations associating long g amma-ray bursts 
(GRBs) with core -collapse supernovae dHiorth et al.l 120031 : 
IStanek et"an l2003|) provid e strong support fo r the collapsar 
scenario ( Wooslev 1993; Paczvnskil [Ull; iPopham et all 
[1999; MacFadyen & Woosley 1999^ In this picfiire, the 
central k,1>Mq of a massive star collapses directly to a black 
hole, while infalling material with higher specific angular 
momentum forms an accretion disk of mass ~1Mq. Viscous 
stresses in the disk drive material toward the hole at rates of 
>O.1M0S"', which powers relativistic outflows responsible 
for both the GRB and the supernova. Although this model 
has many favorable generic features, certain aspects of the 
stellar collapse and disk evolution remain poorly understood 
or largely unexplored. New observational results demand 
augmentation of the simplest versions of the collapsar 
scenario. 

Recent Swift observations have revealed that X-ray flares 
occur ring ftil 0^-10'^ s afte r the initial burst are quite common 
(e.g., O'Brien et al.ll2006l) . The flares exhibit a wide range of 
amplitudes and time scales, and some GRBs show multiple 
eruptions. Explanations for t he flares point to delayed activ- 
ity of the central engine (e.g..lKing et al.P2005|: lBurrows et alj 
20051: iP roga & Zhang 2006; Perna et al. 200l). In particular, 
Perna et al.i 6200 6 ) speculate that an extended collapsar disk 



may be gravitationally unstable and fragment at large radii, 
where the viscous time scale coincides with the onset times 
of the flares. This work prompted us to examine in more de- 
tail the physics of the fragmentation process under the ex- 
treme conditions of collapsar disks. We also take a fresh 
look at the observable consequences, with special emphasis 
on gravitational-wave emission. 

Central to the issue of fragmentation in a self-gravitating 
accretion disk is the cooling rate of the material. If the disk is 
formally unstable, gravitationally bound clumps appear only 
if the cooling time is shorter than t he orbital period about the 
central object (e.g., Gammie 2001); these points are reviewed 
in § 121 Radiation leakage is far too slow in collap sar disks to 
be an effective coolant (e.g.. lNaravan et ani2001l) . However, 
such disks have sufficiently high temperatures (> 1 MeV) and 
densities (>10'^gcm"-') that nuclear processes play critical 
roles in the energy budget. We argue in §[3] that nuclear pho- 



todisintegration is the only energy sink fast enough to permit 
fragmentation. The condition for gravitational instability and 
the strong temperature sensitivity of the photodisintegration 
rate allow us to nicely constrain the location, size, and mass 
of an unstable clump. Our estimates suggest that a fragment 
will collapse until it is supported by neutron degeneracy pres- 
sure, forming a low-mass (« 0.1 - IM©) neutron star. 

As described in § |4l a fragment will be driven toward the 
black hole as disk stresses and gravitational radiation remove 
orbital angular momentum. Migration ends with the tidal dis- 
ruption of the fragment. We calculate the time scales for mi- 
gration and subsequent accretion, and consider whether frag- 
mentation accounts for the X-ray flares in long GRBs. As- 
suming that only a single fragment orbits the central black 
hole, we estimate in §|5]the characteristic gravitational-wave 
strain from this binary using a total migration rate that in- 
cludes viscous torques. 

2. GRAVITATIONAL INSTABILITY IN COLLAPSAR DISKS 

Here we consider the conditions for gravitational instabil- 
ity and fragmentation of ^teaiiy-^rare collap sar disks (see also 
iNaravan et 311120011: iDi Matteo et alj|2002h . where the mass 
inflow rate is constant in radius. The steady-state assump- 
tion is a gross simplification, especially at large disk radii 
dMacFadven & WooslevI 19991) . where material in the collaps- 
ing star with high specific angular momentum may continue 
to rain in for many viscous times. Since this region of the 
disk is the staging ground for our study, our results should be 
viewed as only rough approximations. 

Our background model consists of a central black hole of 
mass Mbh surrounded by a Keplerian accretion disk with 
radially varying orbital frequency il = (GMbh/'"^)'^^, sur- 
face density S, isothermal sound speed Cj, and vertical scale 
height H = Cj /^l. Mass flows at a rate of M = 37ri/I] (e.g., 
iPringle 198 1|), where v= aCsH i s the usual viscosity prescrip- 
tion ( Sh akura & Sunvaevi 11973 ). We assume that H /r = -q 
is a fixed parameter, which is physically plausible in the 
thick, advectiye out er region of a steady collapsar disk (e.g., 
iPopham et al.lll999h . Results below are given in term of the 
scaled variables rioo = ^/(IOOGMbhc"^), = Mbh/3Mo, 
Ml =M/1M0S"', ao.i = a/0.1, and 770.5 

Gravitational instability arises when 



:r j/0.5. 
ffoomrd 119641: 



2 



PIRO & PFAHL 



iGoldreich & Lvnden-Bel]|fT965h 

Q=-4;<Qcn,^l. (1) 
ttGIj 

The value of Qcrit depends on the disk thickness and the 
equation of state. The disk equations above imply that Q oc 
arf'M'^ r"-^/^ decreases with r, and thus instability holds for 

ri^>Aaf,ril,Mf\ (2) 

When Q < Qan, the linearly unstable mode with the fastest 
growth rate («il/0 has a wavelength of KiQH and mass of 
(e.g., § 6.3 of Binney & Tremaine 1987) 

(QHfj: « 0.02e2a-\M3Mir^/>o . (3) 

We identify {QH)^Yi with the mass of a bound clump if cool- 
ing is rapid enough to permit collapse (see §[3]). 

Numerical simulations of thin disks with idealized equa- 
tions of state and cooling prescriptions show that gravitational 
instability leads to fragmentation when the cooling time scale 

satisfies ^ti.„„\ < ^, wher e Q » TO M^'rim^ s~', and $ » 1 - 
10 (e.g., iGamjnie 2001; Meiia et a l. 2005; R ice et al.ll2005h . 
When cooling is slow, gravitational instability produces spiral 
waves that dissipate as thermal energy, driving the disk back 
toward stability. It is unclear if we can reliably use previous 
simulations to draw conclusions about collapsar disks, but we 
take these results as a plausible starting point. 

3. COOLING MECHANISMS 

Using the disk equations in § |2l we evaluate the discrim- 
inant rifcooi = VlTiC^/qH for various possible cooling mech- 
anisms, where q is the relevant cooling rate in units of 
erg s"' cm"^. Cooling via radiative diffusion can be ignored 
outright, since huge optical depths lead to photon diffusion 
times of many years ( lNaravanetal.ll2001h . 

The most important neutrino cooling processes are pair 
annihilation to neutrinos and th e Urea process . We adopt the 
corresponding cooling rates in iPopham et al.l (1 19991) . These 
rates are strong functions of disk temperature, which for an 
ideal gas is T\q = r/10'°K « 3. 6/^,7705-^3 '"loo where /i is 
the mean molecular weight. Here we assume for calculational 
purposes that material at rioo > 1 is pure helium (/i = 4/3) 
and has the ideal-gas equation of state. At the radius where 
Q = Gcrit (see eq. Q), cooling from pair annihilation gives 

(f^fcooi)pa.r ~ SOOaJ/fijo'sM^'M;'^^ , (4) 
A similar calculation for Urea cooling yields 

(rifcool)urca « SOa^' T^o^M^' M^' X^^ , (5) 

where X^uc, the mass fraction of free nucleons, is very sen- 
sitive to T (e.g.. iPopham et al.|[T999l) . Both (fifcooi)pair and 
(JlfcooOurca secm too large to allow fragmentation when ao.i ~ 
1 (a natura l choice for self-gravitating disks; Gammie 2001; 
iRice et akl 12005) and X^uc ~ 1- While we can not rule out 
pair annihilation as a relevant coolant (for instance, when 
ao.i ~ 0.1), Urea cooling is important only where nuclear 
photodisintegration is nearly complete. As we now demon- 
strate, photodisintegration itself will absorb sufficient energy 
quickly enough to allow fragmentation. 
For photodisintegration of "^He we set 

(f^fcooi)photo = ^ « 20r?g.5M3-V-^f A-1 , (6) 



where £ = 7.07 MeV/m^ is the "^He binding energy per nu- 
cleon and A is the reaction rate in s"' . Experiments show that 
the '*He(7,n)^He and "*He(7,p)T cross sections rise abruptly 
at photon en ergies of w20 MeV t o a plateau value of ct w 1 mb 
(10"^'' cm^; iFeldman et al.|[T990l) . Subsequent photodisinte- 
gration steps have larger cross sections and smaller t hresh- 
old energies (ISkopid et al.lll9811 : lBirenbaum et alj 19851) . con- 
tributing little to the overall rate. At temperatures kT < 
20 MeV (Tio < 23.2), we approximate the '^He photodisinte- 
gration cross section as a step function and integrate over the 
Wien tail of the Planck spectrum to estimate 

A«10'^r,oexp(^-^)s-', (7) 

which is a good fit to detailed calculations provided to us by 
H. Schatz (2006, private communication). Because of the ex- 
ponential factor, A varies from «1 to wlO^ as Tio is increased 
from 0.6 to 1.2. Thus, photodisintegration turns on sharply at 
rioo ~ 67705, interior to which (ilfcooi)photo plummets to ^1. 
Rapid photodisintegration may initiate at slightly larger radii 
if a gravitationally unstable clump has some compressional 
heating as it begins to collapse. 

It seems that photodisintegration is a very effective coolant, 
permitting fragmentation over a small range of radii near the 
location where Q = Q„it in the unperturbed disk. By substi- 
tuting the radius at which 2—1 (see eq. ||2l) into eq. (O, we 
estimate the mass of a bound fragment: 

Mf^Q.lTjl^MiMQ. (8) 
Multiple fragm ents may form in the same region (e.g., 
iRice et an i2005). although it is unclear how their individual 
masses are distributed. We assume in § §|4] and |5] that small 
fragments merge into a single body of mass wO.l-lM©. 

The breakdown of "*He into free neutrons and protons 
has two added benefits. First, efficient Urea cooling be- 
comes possible, promoting the continued collapse of the 
fragment. Second, as photodisintegration proceeds and 
A'nuc approaches unity, the adiabatic pressure-density expo- 
nent Fi = (91ogP/91ogp)ad may dip below 4/3, in which 
case the gravitationally bound fragment is unstable and col- 
lapses dynamically. This s econd point is illustrated by 
llmshennik & NadezhinI (1 19661) . who compute the thermody- 
namic properties of a plasma composed of ^^Fe, "^He, neu- 
trons, protons, electrons, and positrons. We have repeated 
these calculations' without ^^Fe. Figure [1] shows the results 
in the temperature-density plane, where we see that the tra- 
jectories of our disk models pass through the region with 
Fi < 4/3. However, the overlap of 2 < 2crit (thick, solid 
curves) and Fi < 4/3 occurs only for sufficiently high M (e.g., 
>lMoS-i). 

Collapse of the fragment halts once the density is high 
enough for neutron degeneracy to provide pressure support. 
At this point we might expect the fragment's radius to be 
Rj ~ 10(M//MQ)"'''-'km. However, there are numerous un- 
certainties regarding the formation and ultimate structure of 
this newly formed low-mass neutron star, including the hy- 
drodynamics of the collapse, the equation of state, the effects 
of the neutrino pressure and finite temperature, and the fate of 
the fragment if Mf < 0.1 Mg, t he minimum mass of a stable 
neutron star (e.g.. lColpi et al.ll989) . An accurate, quantitative 
assessment of these uncertainties is not possible at present. 

' The assumption of thermodyna mic equilibrium is no t strictly valid unless 
the material is opaque to neutri nos IBeloborodo\ll2003h . but the approach of 
llmshennik & NadezhinI (1963) is still useful as an illustrative tool. 
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Fig. 1. — Properties of a plasma composed of ''He, protons, neutrons, 
electrons, and positrons as a function of density and temperature. The thin 
solid lines are contours of constant Fj, and the thin dotted line mai'ks the 
half photodisintegration contour (Xhc = X, -vXp). Thick lines represent disk 
models with a = 0. 1 , 77 = 0.5, and the accretion rates shown. Dashed and solid 
portions of these curves indicate whether the disk is gravitationally stable or 
unstable (Q < Qcrit). respectively. 

4. MIGRATION AND TIDAL DISRUPTION 

A net torque on the orbit of the fragment arises from dissi- 
pation within the disk and the loss of orbital angular momen- 
tum in the form of gravitational waves. When the fragment is 
massive enough to open a gap in the disk, it migrates inward 
on the viscous tirne sca le r/r w r^/iy (i.e., type II migration; 
iLin&Papaloizoull 19861) . or 

« iaTj^ny' « 1 ao,\77o.|M3rJ^o s . (9) 

Gap formation is expected when Mj/M^u > a^^-rf w 

0.08 a^,^ 7^0 5 (e.g.. iTakeuchi et anil996h . a result derived in 
relation to planet formation. For simplicity, we adopt as the 
migration time from disk torques, even if the gap condition 
is not strictly satisfied. Loss of orbital angular momentum in 
gravitational waves causes inspiral on a time (Peters 1 964) 

fgw = ^(^f^) «700Mf/VVfooS, (10) 

where M = M iMq w M^J^M^^ is the chirp mass. The two 
migration times are equal when the orbital frequency is 

3 / «: 2 \ 3/5 

"-"^(^) -4.8aK.^A<-kHz. (11) 

For n > 57eq, gravitational radiation yields the dominant 
torque. Both ilt^ > 1 and fifgw > 1 when r > GMbh/c^, 
so that an initially circular orbit remains nearly circular dur- 
ing migration. Also note that efficient neutrino cooling at 
''ino ^. 1 allows the disk to thin substantially (77 w 0.1-0.3; 
iPopham et al.1 11999'). increasing t^, and decreasing figq. This 
is important in § |5] where we estimate the gravitational- wave 
signal. 



Tidal disruption of the fragment begins when it fills its 
Roche lobe of radius w0.5r(M//MBH)''^^. If the fragment has 
the radius of a neutron-degenerate object (see §[3]) disruption 
occurs at frequency 

r!dis«3.5(My/M0)kHz. (12) 

Whether the fragment is disrupted dynamically or under- 
goes stable mass transfer driven by gravitational radiation 
(e.g., Bildsten & Cutler 1992) , most of the fragment's mass 
is stripped in < 1 s. 

Given the short time for disruption, it seems u nlikely that 
fragmentation of the sort discussed here and in iPerna et aH 
(I2OO6) can account for the late-time X-ray flares in long 
GRBs, which have rise and decay times similar to their onset 
times of >10^s. Moreover, if a fragment forms at rioo ~ 1^ 
and viscous stresses set the migration time to wlO-lOOs, ac- 
cretion of the fragment might blend with the prompt emission. 
The migration timescale may be increased if the fragment is 
espec ially massive in comparison to the disk (Clarke & Syer| 
Il996l) , but probably by a factor of < 10. Better understanding 
of the interaction between the fragment and the disk will re- 
quire detailed hydrodynamical simulations. Gravitational in- 
stability at large disk radii may power the GRB flares, but, if 
so, probably through enhanced angular momentum transport 
rather than fragmentation. 

5. DETECTION OF GRAVITATIONAL WAVES 

LIGO and similar detectors, such as VIRGO, GEO600, 
and TAMA300, are sensitive to gravitational waves with 
frequencies of ~ 10-10^ Hz. For the full range of GRB 
models, gravitational waves in the LIGO band would be 
generated within the central engine, either the compact bi- 
nary progenitor of a short burst or in a collapsar disk. 
Measurement of the gravitational-wave signal from a GRB 
would be a uniq ue probe of the physics at the heart of 
the explosion (se e 'Kochanek & Piran"1993'; ^ Finn et akll 19991 ; 
Ivan PuttenI 12001 ; Kobavashi & Meszaros 20031). Previous 
work on gravitational radiation from collapsar disks have 
addressed gene ral no naxisymmetric instabilities in compact 
accretion tori dvan P utten 2001) or bar- mode instabilities 
dKobavashi & Meszlros .2003) . Unlike these earlier scenar- 
ios, the simplest form of our fragmentation model constrains 
the mass quadrupole moment and its evolution, features we 
now exploit to obtain the gravitational-wave signal. 

As in the previous section, we imagine a binary composed 
of a single fragment in a circular orbit around the black hole. 
The measured gravitational-wave strain is h{t) = h()(t) cos (f>{t), 
where / = 100/iooHz = JI/tt is the slowly varying wave fre- 
quency, and (p = Itt J dtf is the accumulated phase. For a 
source at distance D = lOODiooMpc, the strain amplitude is 

~ 6.4 X 10-24eA^f /f/o'DT^o , (13) 

where 8 = 0-4 is a factor incorporating the orientation of 
the source and the antenna pattern of the detector (e.g., 
iFinn & Chernoffll993 ). For completely random source orien- 
tations and positions on the sky, the root-mean-square of Q is 
■\/ (0^) = 1 .6. Long GRBs seem to be narrowly beamed, prob- 
ably in a direction normal to the accretion disk. Therefore, we 
expect an associated black-hole/fragment binary to have a low 
incUnation, which gives \/(&^ — 2.5, or ~60% larger than 
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Fig. 2. — Characteristic gravitational-wave strain from tlie inspiraling 
black-hole/fragment binary for = 1 00 Hz {thick solid line) and 1 ItHz (thick 
dashed line). Both curves assume M = \Mq, D = 100 Mpc, and © = 2.5. 
The thin lines show the broadband strain noise for current LIGO (dotted), 
enhanced LIGO^ (dot-dashed), and advanced LIGO (dashed). 



the random case (see also iKochanek & Piranlll993l) . In what 
follows, we let 9 = 2.5. 

When the waveform is known, the matched filter approach 
can be used to search for the sig nal. The signal-to-noise for a 
slow binary inspiral is then (e.g.. lFlanagan & HughesllT998h 



I) =/^<.n/)..,(/)|^; 



(14) 



where Wcyc = /^// is the number of cycles per unit In/, and 
hn is the strain noise in a bandwidth centered on /. We see 
that 

'\J ^cyc 1^ 

the signal gain and ^Wcyc ha = he is the associated 
characteristic strain. 

If \r/r\ = tl}+t~l, is the total inward migration rate from §|4l 
then Hcyc = ^fr/'ir. Some minor algebra yields 



5/3- 



H -1 



(15) 



where /eq = f^eq/Tr (see eq. ifTTl ) tunes the importance of vis- 
cous torques. When / <C /eq viscous migration dominates and 

^ For further information see "Enhanced LIGO", LIGO internal docu- 
ment T060156-01-I, R. Adhikari, P. Fritschel, and S. Waldman available at 
[http://www.ligo.caltech.edU/docs/T/T060 156-01 .pdf 



he oc f^/^, while gravitational-wave driven inspiral (/ 3> /eq) 
gives he (X /~'^^. Figure |2] shows heif) for /eq = 100 Hz and 
1 kHz. The lower value of /eq essentially assumes that the disk 
is fairly thin (rj w 0.1), as expected when neutrino cooling is 
efficient. Since he is peaked near /eq, albeit rather broadly, 
we can use the approximation S/N w /j£-(/eq)//^«(/eq)- Inspec- 
tion of Figure|2]indicates that, forDioo ^ 1, 5/A^« 5-10 when 
/eq = 100 Hz, while S/N<1 when /eq = 1 kHz. Detection with 
advanced LIGO seems quite possible for nearby GRBs. The 
combined effects of viscous torques and gravitational-wave 
emission on the orbital evolution create a distinct strain sig- 
nal, where an estimate of /eq would provide unique insight 
into the physics of collapsar disks. 

Estimates of the local rate of detectable long GRBs span 
wO.1-1 Gpc"^ yr-i (e.g.. lGuetta et alJl2005h . based largely on 
a small sample with moderate to high redshifts. This implies 
a rate of <4 x lO'-'DjgQyr"' within a distance D . However, 
two GRBs in the past 8 years within 150 Mpc (980425 and 
060218; e.g., Galamaetal. 1998; Ferrero et al. 2006) sug- 
gest tentatively that the local rate is more uncertain and could 
be as large «0.1 DjQQyr"', which is consist ent with 10"* of 
the lo cal rate of core-collapse supernovae (ICappellaro et alJ 
Il999l) . Beaming implies there may be wlOO unseen GRBs 
for every detectable event, giving a total local GRB rate of 
<10Z)j()Qyr"'. If a misdirected GRB is detected by multiple 
gravitational-wave antennae, it may be possible to locate the 
source to within a few degrees and provide a target region to 
search for an "orphaned" electromagnetic afterglow. 
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